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ABSTRACT 

In a concordant A Cold Dark Matter (ACDM) model, large-angle Cosmic Microwave 
Background (CMB) temperature anisotropy due to linear perturbations in the local 
universe is not negligible. We explore a possible role of an underdense region (void) 
that may cause an anomalous Cold Spot (CS) in the CMB sky. Although the observed 
anomalous cold region with a surrounding hot ring can be produced by an underdense 
region surrounded by a massive wall, a decrement in the CMB temperature in the 
line-of-sight is suppressed because of blueshift of CMB photons that pass the wall. 
Therefore, undercompensated models give better agreement with the observed data 
in comparison with overcompensated or compensated models. We find that it is likely 
that ~90 per cent of the CMB fluctuation is generated due to an overdense region 
surrounded by an underdense region at the last scattering surface, and the remaining 
~10 per cent is produced due to a single spherical underdense region with a radius 
r ~ 6 X 10^ h^^Mpc and a density contrast 6m ~ —0.009 {2a) at redshift z 1 in the 
line-of-sight to the CS. The probability of having such two aligned structures is ^ 0.7 
per cent if the perturbed region at z ~ 1 is moderately undercompensated. 

Key words: cosmic microwave background - cosmology - large scale structure of 
the universe. 



1 INTRODUCTION 

One of the interesting non-Gaussian features found in the 
cosmic microwave background (CMB) sky is the p resence of 
the C old Spot (CS) in the Southern hemisphere (jCruz et al.l 
l2005h . Although the statistic al significance is not signifi- 
cantly high (1.85 per cent) (|Cruz et al.l |2007| ). it maybe 
a hint of new physics that is absent in the concordant 
ACDM model. The origin of the temperature decrement 
may be a time-varying gravitational potential of a super- 
void wit h radius 200 — 300/t~^M pc in the line-of-sight to 
the CS (|lnoue fc Silkll2006l . 120071 ') . or t hat of a cosmolog - 
ical texture at high redshifts z > 1 (|Cruz et al.l |2007| ). 
So far, no clear ev idence of supervoids in the line-of-sight 
has been obta i ned jRudnick et al]|2007l : I Bremer et al. I I2OIOI : 
ICranett et al] l201oi r" If such a huge supervoid exists, it 
will affect various kinds of observables (|Das fc Sperge]|[2009l : 
iMasina fc NotarillioOQbl lal'). Therefore, it can be easily dis- 
tinguished from other scenarios such as a cosmic texture 
l|Cruz et al.ll200'7l ) in the near future. 

A key issue is an evaluation of the effect of non- 
compensation of mass. It has been pointed out that the CS 
is anomalous because of a hot ring-like str ucture around it 
IZhang fc Hutereijl20ldl : llnoue et al] I2OI0I '). However, pro- 
ducing a hot ring by a linear or quasi-linear compensated 
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void is difficult due to the linear integrated Sachs-Wolfe 
(ISW) effect though such a feature c an be realised by non- 
linear compensated c lusters or voids (iTomita fc Inouell2008l : 
ISakai fc Inou3 12OO8I : llnoue et all |201GI : IZhanll201ll ). If one 

considers a massive wall that overcompensates the low den- 
sity region inside, such a feature might be easily made even 
in the case of linear fluctuations. However, we need to esti- 
mate the chance of having such structures. 

Moreover, in order to explain the entire feature of the 
CS by a single supervoid, one has to assume an exception- 
ally huge void, which is dif ficult to produce i n the stan- 
dard cosmological scenarios (jlnoue fc Silkll2007l ). In fact, it 
is more natural to assume that most of the feature of the 
CS is produced by fiuctuation s in the gravitation al potential 
at the last scattering surface (|Valkenburdl2012l ). Therefore, 
one must take into account various effects such as the ordi- 
nary Sachs- Wolfe (OSW) effect (due to spatial fluctuations 
of the gravitational potential and the temperature) and the 
Doppler effect due to the bulk velocity of photon-baryon 
fluid at the last scattering epoch. It is of crucial importance 
to study the expected value of the percentage of the contri- 
bution from an assumed supervoid (or a low density region) 
to the CS and the most probable redshift in the framework 
of standard scenarios that predict Gaussian primordial per- 
turbations. Even if the CS has been produced as a fluke, an 
estimation of a possible contribution from the local struc- 
tures in the line-of-sight to the CS is important since the 
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ISW contribution is expected to be maximum in that direc- 
tion. 

The present paper is organised as follows. The linear 
and the non-linear ISW effects due to an non-compensated 
void (or cluster) are explored in section 2. The statis- 
tical significance of a moderately undercompensated un- 
derdense region that will partially explain the CS is dis- 
cussed in section 3. We conclude the paper in section 
4. In the following, unless noted, we assume a concor- 
dant ACDM cosmology with {Qrn,,o,^A,o,^b,o,h,as,n) — 
(0.26,0.74,0.044,0.72,0.80,0.96), which agre es with the re- 
cent CMB and large-scale structure data l| Sanchez et al.l 
l2009l ). 
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2 MODEL OF NON-COMPENSATED 
VOID/CLUSTER 

In order to investigate temperature variation due to non- 
compensated underdense or overdense regions, we consider 
spherically symmetric density fluctuations with a comoving 
radius ro surrounded by a wall at ro < r ^5 ri. To localise 
the perturbation, we require that the perturbation is glob- 
ally compensated. To do so, we introduce a cutoff comoving 
radius r2 ^ ri where the total mass of the perturbation is 
compensated if integrated over < r < r^. In what fol- 
lows, for simplicity, we only consider top-hat type matter 
density perturbations for void (cluster) whose linear density 
contrast is given by 

( (5o, r<ro 
S"^ = I St, ro^r<ri (1) 

where 5q , 5i and $2 are constant (Fig. 1). From a condition 
that the perturbation is local, we have Sgrg + di{r'l — rg) -f 
^2'(f"2 — rf) = 0. To describe the local mass compensation of 
the perturbation, we define a parameter e as 

_ ^o^rg + ^f(r!^rg) 

e = corresponds to a locally compensated case, whereas 
e > 0, (or e < 0) corresponds to a locally overcompensated 
(undercompensated) case in the region defined by ^ r ri 
in linear order. Even in the case e = 0, if ri — ro > ro and 
n = r2, we call it a " moderately over(under)compensated" 
model in the following]. 

In order to incorporate non-linearity of the scalar- 
type perturbation, we use second order perturbation theory 
i Tomita fc I nouc 2008). In what follows, we assume that the 
cosmological Newtonian approximation holds for the pertur- 
bation that we consider. In other words, we assume that the 
peculiar velocity is sufficiently smaller than the light veloc- 
ity, and the size is sufficiently smaller than the present hori- 
zon scale. In what follows we only consider growing mode 
of density perturbation that is adiabatic. The temperature 
anisotropy due to time-evolving gravitational potential per- 
turbation ^ is given by an integration along with a geodesic 
of the CMB photon that passes through the potential. 




Alternatively, we may call it "a model with a thick wall" . 



Figure 1. Linear matter density contrast for a top-hat type 
spherical underdense region. 

where rj is the conformal time. If the amplitude of the den- 
sity contrast S is sufficiently small, the gravitational poten- 
tial can be expanded as a sum of linear order and second 
order perturbations ip = ip^ + . Time evolution of the 
linear perturbation tp is written as a product of a potential 
function F of comoving coordinates x and a growing factor 
D{ri), which is a function of the scale factor a(ri) as 

V'^(r,,x) = -i3(r?)F(x) 

= i(l-^P')i.(x), (4) 

where 

oilmO 

X / daa^^^in^o+^Aoa]-"-^^ + -^a. (5) 

Jo imO 

Here, fimcf^Ao are density parameters of non-relativistic 
matter and a cosmological constant A of the background 
FRW universe and a prime denotes derivative with respect 
to the conformal time rj. In a similar manner, the second 
order perturbation can be written as 

^^(»7,x) = UV)FA^)FA^) + 1006(r?) • *o(x), (6) 
where is dF/dx\ and 

and ^'o and Q{ri) satisfy 

Q" + ^Q' = -[p-liP')% (10) 

respectively (Tomita & Inoue 2008). First, we consider ef- 
fects of linear perturbations. As shown in Fig. 2, a hot ring 
that surrounds a cold spot appears as the wall around an 
underdense region becomes more massive than a mass defi- 
ciency inside since a decay of gravitational potential at an 
accelerated epoch leads to an increment in the temperature 
of the CMB photons. Similarly, a cold ring that surrounds 
a hot spot is produced as an underdense region surround- 
ing an overdense one becomes less massive than a mass that 
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Figure 2. Effect of noncompensation: temperature anisotropics 
as a function of angular radius d from the center of the pertur- 
bation at redshift z = 1 are plotted for voids (left) and clus- 
ters (right). The parameters of the mass distribution are set to 
(ro,ri,r2) = {0.1H~^ ,Q.lbH~^ ,O.ZH~^). Note that the vertical 
scales of the upper and lower panels are different. 



Figure 4. Effect of noncompensation; profiles of gravitational po- 
tential for non-linear density perturbations at 2 = 1 are plotted 
as thick curves and those for corresponding linear density per- 
turbations are plotted as thin curves (z = 1) and dashed curves 
(z = 10) as a function of comoving distance r from the center. 
The parameters are the same as in figure 3. 
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Figure 3. Effect of noncompensation: profiles of mass density 
contrasts for non-linear density perturbations at z = 1 are plot- 
ted as thick curves and those for corresponding linear density per- 
turbations are plotted as thin curves (z = 1) and dashed curves 
(z = 10) as a function of comoving distance r from the cen- 
ter. The unit of the length is the present Hubble radius . 
The linear density contrasts at the center at z = 1 are as- 
sumed to be either (5g = —0.6 (left panels) or 5^ = 0.6 (right 
panels). The mass is either locally compensated (upper panels), 
overcompensated (lower left) or undercompensated (lower right). 
The parameters of the mass distribution are set to (t-q ,r\,r2) = 
(0. \ 0. 1 Si^o" \ 0. 3Ho" ^ ) . 



point r = r\, the gravitational potential becomes concave 
outward. For this reason, the neighbourhood of r = ri con- 
tributes a negative mass to the density distribution, which 
leads to a shift of the wall towards the outer (inner) shell 
r > r2 in the case of a void (cluster). On the other hand, an 
additional mass due to the second order efTect leads to a de- 
crease in the gravitational potential in the inner low density 
region ^ r < ri. Thus, the gravitational potential expands 
(shrinks) faster than the corresponding linear potential in a 
void (cluster) in the outer region, but the value of the poten- 
tial decreases in the inner region as shown in Fig.43. These 
changes lead to a temperature decrement (increment) in the 
line-of-sight to the neighbourhood of the center (edge) of 
the perturbation regardless of the sign of the fluctuation as 
shown in figure 2. 

For non-compensated quasi-linear perturbations, both 
of the features appear. As shown in figure 2, for overcom- 
pensated quasi-linear voids, we would observe a hot ring that 
surrounds a prominent cold spot. In contrast, for undercom- 
pensated quasi-linear clusters, we would see a cold ring that 
surrounds a cold spot. 

Thus it is important to study locally overcompensated 
(linear/quasi- linear) voids in order to explain the origin 
of the Cold Spot because it seems to have a hot ring 
dlnoue et al.ll201(ll ). 



compensates the overdense region inside. These behaviors 
agree with the predicted temperature variation due to the 
linear ISW effect. 

Next, we consider quasi-linear regime. Although the 
non-linear mass seems to be overcompensated at r = r2 
with £ = (see Fig. 3), it is actually locally compensated. 
As shown in Fig.4, the spatial derivative of the potential is 
not continuous at ri. Therefore, it gives a contribution that 
is proportional to Dirac's delta function Soir — r\). Fur- 
thermore, if it is smoothed out in the neighbourhood of a 



3 STATISTICAL SIGNIFICANCE 

In order to estimate the statistical significance of our model, 
we use CMB temperature anisotropy AT/ that is filtered by 



^ Interestingly, for quasi-linear clusters, one can see that cancel- 
lation between the linear and second order ISW effects can occur 
and the net temperature anisotropy can vanish near the center if 
the perturbation is compensated at linear order. 
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a spherical top-hat type compensating filter defined as 



{6 < e,„) 

'1 {din ^ d 6 out); 



(11) 



where 9out = cos~^ (2cos&i„ — 1). 

In the harmonic space, the variance of filtered temper- 
ature AT/ can be written in terms of the angular power 
spectrum Ci as 

I 

where A{9in) = 27r(l — cos6li„) and 



w, = 



+ 1) 



/ + 1 



+ XoutPl{Xout) - Pl-l(Xout) 



(13) 



where Pi is the Legendre func tion. 

For detail, see appendix of llnoue et al.l(|2010h . Using the 
WMAP 7-year data, a deviation of Ar/(CS) = -72fiK is 
found at smoothi ng scale 9in = 12° if centered at the center 
of the Cold Spot (jl noue et al]|2010l 'l. This corresponds to a 
posteriori significance level of 4a in the concordant ACDM 
model, since la{9i„ — 12°) = ISfiK for a single realisation. 

Firstly, we assume that a presence of an underdense 
region in line-of-sight to the Cold Spot contributes 18/iK 
deviation (la) of the filtered temperature at smoothing 
scale 9in = 12° as an ansatz. The effects of temperature 
anisotropy due to fluctuations at the last scattering surface 
or its neighbourhood such as the ordinary Sachs- Wolfe ef- 
fect, the early integrated Sachs- Wolfe effect and the Doppler 
effect will be considered in next section. 

In order to assess the effect of non-compensation of 
the surrounding mass, we need to flx shape parameters 
(ro,ri,r2), degree of overcompensation e, and redshift z of 
the center. In terms of comoving distance dp{z) to the center 
of the perturbation, the subtending angles of shape param- 
eters are given by 9i = arctan(ri/cip(2:)), i = 0, 1, 2. For sim- 
plicity, we consider four types of models by flxing the outer 
boundary r2 or 82, the next inner boundary ri and e as 
shown in table 1. Then we calculate the linear density con- 
trast 5q{z) at the center that will give a filtered temperature 
anisotropy 18/xK as a function of the inner most boundary 
ro and redshift z. In order to assess the probability of a 
single realisation of such a linear perturbation, we divide 
the obtained amplitude of the linear density contrast by the 
standard deviation o-f{z) of a top-hat type matter density 
perturbation expected in the concordant ACDM model we 
defin ed. We use an analyt ic form of the CDM transfer func- 
tion (|Bardeen et al.|[l98 6l) and that of the growth factor for 
the Friedmann-Robertson- Walker background with dust and 
a cosmological constant. 

As shown in Fig. 5, these 4 types of models give sim- 
ilar values for the best fitting parameters 80 = 14° — 15° 
and z ~ 1.0. Provided that the perturbation is at linear 
regime, the temperature ani sotropy due to the linear ISW 
effect is proportional to 5o rp l|lnoue fc Silkl2006l ). Therefore, 
objects at higher redshifts will give larger amplitudes in the 
temperature anisotropy for a fixed density contrast and a 
subtending angle in the accelerated epoch. However, if the 
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Figure 5. Significance of realisation of a single top-hat type un- 
derdense perturbation that will give a signal ATj- = —18fiK{la) 
for type A(upper left), type B(uppcr right), type C(lower left), 
and type D(lower right). Contour maps of the amplitude of the 
linear density contrasts Sq at the center in unit of the standard 
deviation ufiz) in a concordant ACDM model are plotted as a 
function of subtending angle of the inner boundary 80 and redshift 
z of the perturbation. 



redshift is high enough to decelerate the cosmic expansion, 
then the linear ISW effect is suppressed. Thus the presence 
of peak at 2: ~ 1 can be explained as a result of these two 
different effects. As long as the perturbation is at the linear 
regime, the place of the peak does not change because the 
temperature anisotropy is proportional to a density contrast 
at the center ^o". 

Although, a hot ring can be produced by a massive wall 
(e > 0) that overcompensates the underdense region inside, 
overcompensated models (type B and C) are less probable 
in comparison with a compensated model (type A). For in- 
stance, the minimum values (most probable) for type A, 
B, and C are 5.6af ,6.1aj , and 6.8a^ . For overcompensated 
models, as we have seen in the previous section, the height of 
the gravitational potential at the center decreases if the sign 
of the potential at the wall is negative (Fig. 4) provided that 
the linear density contrast at the center is fixed. Therefore, 
redshift effect near the center becomes much weaker and 
blueshift effect at the surrounding massive wall suppresses 
further the overall redshifts of the CMB photons that pass 
through the perturbed region. Thus the linear ISW effect 
for overcompensated models is suppressed. In contrast, the 
height of the gravitational potentials for undercompensated 
models or moderately undercompensated models (such as 
type D) increases in comparison with the other compensated 
models with a thin wall (such as type A) since the height 
of the potential at the center increases. Therefore, the ISW 
effect is enhanced for undercompensated models. 

Second, in order to estimate the statistical signifi- 
cance of our single void (underdense region) model in the 
concordant ACDM model, we need to estimate the CMB 
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Table 1. Parameters of top-hat type underdense region 



Type 


e 


ri 






category 


A 





(fo - 


l-r2)/2 


30° 


compensated 


B 


0.5 


(fo - 


l-r2)/2 


30° 


overcompensated 


C 


1 


(fo - 




30° 


overcompensated 


D 





r2 




30° 


moderately undercompensated 



anisotropy due to other effects such as the ordinary Sachs- 
Wolfe effect due to spatial fluctuation of gravitational po- 
tentials and temperature at the last scattering. As we have 
seen in previous section, the most probable redshift of the 
underdense region is 2; ~ 1. Therefore, we divide contri- 
bution from the ISW effect into two categories, 'early' and 
'late' ones. In what follows, we call the ISW effect due to 
fluctuations of gravitational potential at z > 1.5 the 'early' 
type and at z < 1.5 the 'late' type, respectively. The origin 
of the early type ISW effect is a decay of potential due to 
residual radiation at the time of the last scattering. 

For simplicity, we also assume that the late ISW ef- 
fect mainly comes from a single void (underdense region) 
in the line-of-sight to the CS. In this calculation, the nor- 
malization of the angular power Ci is set by erg = 0.80. 
In order to calculate the late type ISW effect, we numeri- 
cally solve a second-order differential equation for the New- 
tonian p otential $ on la rge angular scales. We also use CMB- 
FAST (Seliak fc Zaldarriaga 1996 ) for calculating the overall 
small-angle power in which the effect of acoustic oscillation 
is non- negligible. However, we neglect the effect of reionisa- 
tion since we are mainly interested in large angular power. 
Once the angular power C; of the CMB anisotropy is ob- 
tained, the standard deviation a = y^jATyp) of the filtered 
temperature anisotropy can be calculated straightforwardly 
using equation (12). As shown in Fig. 6, the overall ampli- 
tude (rms) of the filtered temperature is AT/ = 18 — 25/iK 
for 6in = 4° — 16° . Because of anti-correlation between fluc- 
tuations due to the OSW effect and those due to the early 
type ISW effect, the overall amplitude at large angular scales 
is smaller than the sum of the two kinds of amplitudes. In 
contrast, correlation between fluctuations due to the OSW 
effect and those due to the late type ISW effect is small if the 
corresponding size of the fluctuations are sufficiently smaller 
than the present horizon scale. Contribution from the late 
type ISW effect to the overall amplitude is about 25 per cent 
for 8in — 12° . On smaller angular scales, the contribution be- 
comes smaller. For instance, we flnd that contribution from 
fluctuations at redshift 0.45 < z < 0.75 is just ~ IfiK. This 
result is consistent with the previous calculation for su per- 
structures in the SDSS-LRG catalog (|lnoue et al.ll2010D . 

Finally, in order to calculate the statistical significance 
of a single void (underdense region) in the line-of-sight to 
the CS, we assume that it is described by a moderately 
undercompensated top-hat type linear perturbation with 
01 = 62 ~ 30 degree (type D). We calculate the inner- 
most radius of the perturbation ro that is most proba- 
ble for a given redshift z of the center for a given ampli- 
tude of the filtered temperature AT/ (void) for 9in = 12° 
as shown in the previous section. Then for each spherical 
shell with comoving radius dp{z) — ro < r < dp{z) + ro, 
the probability Pvoid(ATf (void); z) of finding such a fluc- 
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Figure 6. Root-mean-square values of the amplitude of the fil- 
tered CMB temperature anisotropy in a concordant ACDM model 
due to the Doppler effects (blue thick curve), the OSW effect 
(purple thin curve), the ISW effect (green thin curve), the early 
type ISW effect (green dashed curve), the late type ISW ef- 
fect (blue dashed curve), and the ISW effect from a region at 
0.45 < z < 0.75 (purple dot-dashed curve). 
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Figure 7. Overall probabilities as a function of contribution from 
a low density region (void) for various redshifts. 



tuation is calculated based on linear perturbation theory 
assuming the CDM power spectrum (Bardeen et al. 1986). 
We approximate the probability distribution of the filtered 
temperature by a Gaussian distribution since the filtered 
temperature is described as a sum of a number of inde- 
pendent modes in the harmonic space. The number of in- 
dependent perturbed low density region for a given angu- 
lar scale ro/dp{z) is approximated as A*' ~ dp{z)^ /{2rg). 
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Then the probabihty distribution function of the filtered 
temperature fluctuation due to a single void is Pvoid = 
1 — F{x)'^ where F{x) is the cumulative distribution for 
the Gaussian distribution g{x) — 2/^/2% exp{—x^ /2) and 
X = ATf (void) l^{{ATf (void) Y). The probability distribu- 
tion function Plss of finding filtered temperature anisotropy 
|AT/(LSS)| = (72-|Ar/(void|))/iK for = 12° due to the 
OSW, early type ISW and Doppler effects can be calculated 
as Plss = g{y), where y = AT^(LSS)//{(AT>(LSS)F)- 

Then the overall probability distribution function of 
finding such configuration is given by Pan = PLSsPvoid- Note 
that we have neglected any contributions due to fluctuations 
at redshifts z < 1.5 except for the assumed low density re- 
gion (void). This is verified if the center of the low density 
region is at z ~ 1 because the late type ISW effect is maxi- 
mum at this redshift for the assumed angular scale Oq ~ 14° 
as we have seen in the previous section. 

In Fig.7, the overall probability Paii(^ |Ar/(CS)l = 
72/iK) as a function of percentage of contribution from an 
assumed low density region (void), AT/ (void)/72/iK x 100, 
is plotted for various redshifts of the region. We find that 
the most probable values of redshift and contribution to the 
filtered temperature from a low density region are z = 1.0 
and 11 per cent, respectively. The remaining 89 per cent 
of temperature anisotropy is produced by an overdense and 
surrounding underdense region at the last scattering surface 
in the line-of-sight to the CS. The maximum value of the 
overall probability of our model is Pan = 0.68 per cent. This 
result indicates that a contribution from the late type ISW 
effect due to an underdense region in the line-of-sight to the 
CS is not completely negligible. Indeed, it turns out that 
the possibility of having such a contribution is maximally 
7 times larger in comparison with the case without a single 
underdense region in the line-of-sight to the CS. 

The most probable parameters of our model is summa- 
rized in table |3] The temperature profile due to the most 
probable underdense region is plotted in Fig. 8. The ex- 
pected temperature decrement in the direction to the center 
is — 21/jK. The density contrast is Sq = —0.0085 correspond- 
ing to 2.1(7 or a probability of 1 in 10 for a single realisation 
of negative density fluctuation ai z = 1. Therefore, it is not 
rare to have such an underdense region. As we can see from 
Fig. 8, a hot-ring like structure around the cold region can- 
not be attributed to the low density region at z ~ 1. If it is 
attributed to the OSW effect, we need to assume a potential 
hill at the last scattering surface. Therefore, it is most likely 
that a low density region surrounding a high density region 
resides at the last scattering surface. 



4 CONCLUSION 

We have explored a possible role of an underdense region 
(void) that may cause an anomalous Cold Spot (CS) in the 
Cosmic Microwave Background (CMB) sky located at the 
southern Galactic hemisphere. Although we have found that 
the observed cold spot with a hot ring can be reconstructed 
by an underdense region surrounded by a massive wall that 
overcompensates the mass deficiency inside, the CMB tem- 
perature decrement in the line-of-sight to the center due to 
the linear integrated Sachs- Wolfe (ISW) effect is suppressed 
because of blueshift of CMB photons at the wall surrounding 
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Figure 8. Temperature anisotropy due to a decaying potential of 
a single low density region (model D) (red dashed curve) and the 
observed Cold Sp ot (WMAP 7) smoothed i n azim uthal direction 
(blue full curve). For detail, see llnoue et al.ll|201(]h . The horizontal 
line represents the angular radius from the center of the CS. The 
parameters of the low density region are given in table 2. 

an underdense region. Therefore, undercompensated models 
give better agreement with the observed data in comparison 
with overcompensated or compensated models. It is likely 
that ~90 per cent of the CMB fluctuation is generated at 
an overdense region surrounded by an underdense region at 
the last scattering surface, and the remaining ~10 per cent 
is produced due to a single underdense region that is mod- 
erately undercompensated with radius r ~ 6 x 10^ /i~^Mpc 
and a density contrast Sm ~ —0.009 (corresponding to ~ 2a 
for a single realisation) at redshift z ~ 1. The probability of 
having such aligned structures is ~ 0.7 per cent if the un- 
derdense region is moderately undercompensated. The tem- 
perature decrement due to the ISW effect is expected to 
be ~ — 20/iK in the line-of-sight to the center. It should be 
emphasised that the chance of having such an underdense 
region is 7 times larger in comparison with the case with- 
out any underdense regions in the line-of-sight to the CS 
in the concordant ACDM model with Gaussian primordial 
perturbations. If such a linear underdense region does not 
reside at redshift z ~ 1, then either an abnormally large 
quasi-linear void at low redshift z < 0.3, or non-Gaussian 
density fluctuations at z ^ 1 is necessary for explaining the 
origin of the CS. Therefore, it is very important to probe 
the density profile in the line-of-sight to the CS at redshifts 
0.3 < 2; < 2 subtending an angular radius 8 > 15° in f uture 
surve y. Ongoing projec ts such as LSST (jlvezid lioOSl ) and 
SKA (jBlake et al.ll2004 l might be helpful for this purpose. 

We have not considered effects of non-sphericity of the 
density profile and chance alignments between other large 
scale structures in the line-of-sight, which are considered 
to be minor effects. However, it will certainly improve the 
accuracy of model prediction if we can take into account 
these effects as well. These effects will be investigated in our 
future work. 
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Table 2. The most probable parameters for the type D underdense region 



z 5^ 9o 


9i 


ro ATy(void) 


Pall(%) 


1.0 -0.0085 14° 


30° 


570/i-iMpc 7.7/iK 


0.68 
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